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ABSTRACT 


Local heat transfer distributions in impingement cooling have been 
measured by Kreatsoulas [1] and Preiser [2] for a range of conditions which 
model those in actual turbine blades. Including the effects of rotation. 
These data were reported as local Nusselt niunbers, but referred to coolant 
supply conditions. By means of a channel flow modeling of the flow in the 
supply and impingement passages, the same data are here presented in terms 
of local Nusselt number distributions such as are used in design. The 
results in this form are compared to the nonrotating impingement results of 
Chupp [3] and to the rotating but non-impingement results of Morris [4]. 
Rotation reduces the mean Nusselt numbers from these found by Chupp by 
about 30 percent, and Introduces Important radial variations which are 
sensitive to rotation and to leading edge stagger angle. 
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NOMENCLATURE 


English Symbols 


A 

cross-sectional area 

(sq. m) 

Cp 

specific heat under constant 

(J/Kg K) 

d 

Impingement jet diameter 

(m) 

D 

leading edge diameter 

(m) 

de 

equivalent diameter (A/s) 

(m) 

f 

friction factor 


h 

heat transfer coefficient 

W/sq.m K) 

k 

thermal conductivity 

(W/m K) 

M 

Mach number 


Nu 

Nusselt number 


P 

pressure 

(N/sq.m) 

P 

jet hole pitch 

(m) 

r 

position vector 


r 

radius 

(m) 

Re 

Reynolds number 


S ' 

perimeter 

(m) 

S 

surface area 

(sq.m) 

S 

jet span 

(m) 

St 

Stanton number 


Ts 

average surface temperature ^ 

(K) 

Tf 

average fluid temperature 

(K) 

Tt 

total fluid temperature 

(K) 
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Ty, effective wall temperature 
u velocity vector 

u velocity component in r direction 
Vj impingement jet velocity 

z wall to jet distance 

Greek Symbols 

y Cp/Cv = specific heat ratio 

y stagger angle 

viscosity 
p density 

SI angular velocity 

Subscripts 

1 impingement channel 


(K) 

(m/sec) 

(m/sec) 

(ra) 


(rad) 

(Nsec/sq.m) 

(kg/cu.m) 

(rad/sec) 


2 


supply channel 
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I . INTRODUCTION 

One of the most challenging tasks for advancing gas turbine technology 
is to improve gas turbine efficiency by Increased turbine inlet temperature, 
while maintaining low enough metal temperature levels to achieve acceptable 
component life. To maintain acceptable temperatures of turbine components, 
it is necessary to develop effective cooling methods, and the Importance of 
using coolant air in airfoils has been repeatedly emphasized. One of the 
most effective cooling methods for local internal cooling of turbine blades 
is jet Impingement, which is achieved by blowing cooling air through a 
series of holes located Inside the blades. The effectiveness of this 
internal cooling technique is measured by the heat transfer coefficient 
which gives the heat flux per temperature difference between the airfoil 
and the coolant. 

The accurate prediction of the local convection heat transfer coeffi- 
cient is essential in order to optimize the use of the coolant supply and 
minimize thermal stresses. Although numerous Impingement heat transfer 
studies have been conducted and reported, most of the available data have 
been measured in stationary systems, so they do not address the rotational 
effects in impingement cooling systems. Recently, an extensive set of data 
on impingement cooling with rotational effects has been obtained by Kreat- 
soulas [1]. These experiments have revealed large effects due to rotation, 
on both the mean Nusselt number (averaged on the leading edge surface) and 
on the detailed heat transfer distribution, which is created by the imping- 
ing jets. All of Kreatsoulas' experiments were conducted with a -30 degree 
stagger angle, typical of the inlet of a high work rotor at the mean radius. 



10 


modeled as channel flows, using friction coefficients and heat transfer 
coefficients which represent averages around the circumference of the 
passages, and are functions of the radial coordinate along each passage. 

Because of the complexity of the flow. It has been possible to describe 
the effectiveness of Impingement cooling only with reference to specific 
geometries. Thus, Ref. [3] gives correlations In terms of the jet Reynolds 
number, jet diameter, and the geometric parameters. We describe the 
effects of rotation relative to this correlation. 

Some experiments have been conducted in simple rotating cylindrical 
passages which show the effects of rotation In such simple geometries and, 
again, we have compared the data of Kreatsoulas and Preiser with this data, 
to show the effect of the Impingement jets on the heat transfer In the 
rotating passage. Within the channel flow approximation adopted here, the 
effects of rotation which are encompassed in the empirical correlations, at 
least In first approximation, are: 1) the radial pressure gradient in the 
supply and Impingement passages, and its dependence on the heat addition in 
each passage, which is important in that it controls the flow rates in the 
impingement orifices, 2) the effect of buoyancy on the heat transfer, and 
3) the averaged effect of Coriolis induced cross flow in the passages. 
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II. EXPERIMENTAL APPARATUS 

The experimental apparatus for study of Impingement cooling in a 
rotating system was mainly designed and developed at the MIT Gas Turbine 
Laboratory and is described in detail by Kreatsoulas [ 1 ] . A brief descrip- 
tion of the apparatus will be discussed in this section in order to present 
the modeling of the flow in the impingement region. 

The apparatus consists of the rotor, the supporting structure, heat 
exchangers, instrumentation, calibrating body and blade model, as illustrated 
in Figure [1]. The blade model, which simulates the leading edge geometry 
of a real turbine blade, rotates in a vacuum chamber. The vacuum chamber is 
used to reduce heat transfer by convection from the outer wall to the 
environment, as well as to eliminate any absorbing medium between radiometer 
and measurement spot. The Infrared radiometer measures the external 
surface temperature of the foil with a spatial resolution of about 1mm. 

The local heat transfer coefficients throughout the leading edge region of 
the blade are then determined from the measured electrical input power and 
the skin temperature distribution. 

Figure [2] shows a geometric description of the flow system of the 
test model. A thin stainless steel foil (Kanthal A-1), a high resistivity 
and low thermal coefficient of resistivity material, is electrically heated 
by dissipation in the skin while being cooled from one side by impinging 
jets of refrigerant gas (Freon 12). The coolant enters into the supply 
pleniim, flows through the jet holes to the impingement section, and 
exhausts radially to the exhaust channel. Inside the cooling passage of 
the model, type E thermocouples and pressure probes are placed to measure 
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III. FORMULATION OF THE MODEL 

Figure [3] Illustrates a tube rotating about an axis with a constant 
angular velocity. Fluid flows in the tube and the motion of this fluid is 
referred to as a reference coordinate system (r,e,z). The acceleration 
vector of a fluid particle, which has a position vector r and a velocity 
vector u, relative to the rotating coordinate system is 

-► Du 

a = — + 2£lxu + nx(flxr) (1 ) 

-¥■ ->■ ->■ “>■ 

Du/Dt, 2dxu, and S)x(Dxr) refer respectively to the total derivative of the 
velocity vector, Coriolis acceleration, and the centrifugal acceleration. 
The Navler-Stokes equation may then be approximated as 

Du 2-+- ... 1 , 

p TTr + 2nxu + flx(nxr) = -Vp + pV u (2) 


where p, p, and p are the density, viscosity, and pressure of the fluid 
respectively. For a steady flow in the radial direction, Eq. (2) may be 
simplified as 


3u 2 

pu pD r 

3r 


3r 


2 

pV U 


(3) 


Equation (3), however, does not Include the effects of the Incoming 
jets from the supply channel, which we shall assume enter the impingement 
passage with zero radial momentum. If an elementary radial volume, as 
shown in Fig. 4, is considered, the momentum balance becomes 


_ ^ 

dr 


du d , , 

P" to ^ " d7 ‘P"’ 


(4) 


From continuity. 


dA, 


to ‘P"^^ = P'^j dT 


(5) 
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The energy equation canr be derived'ln a similar manner as the equation 
of motion by taking a radial element. The heat transfer rate across the 
surface dS is written in terms of the heat transfer coefficient and the 
temperature difference as 

dq = h(T - T)dS (14) 

w 


As a first approximation, the heat transfer rate is estimated by 
taking the temperature difference between the average temperature of the 
fluid T and the effective wall temperature T^, which is assumed to be 
constant along the passage. T„ is defined as 


T 

w 


^^s " Vf 


(15) 


where 

A<| impingement channel cross-sectional area 
A 2 ° supply channel cross-sectional area 
Tg s average surface temperature 
Tf = average fluid temperature 

In the final heat transfer coefficient calculation, T^ is replaced 
with the measured local wall temperature, and T is replaced with the local 
fluid temperature. This is to be described in detail in IV. 

The average heat transfer coefficient, as a first approximation, is 
estimated by applying the following analogies. The Reynolds analogy, which 
assumes complete similarity of momentum and heat transfer, provides the 
following equation. 

Nu ^ h 

Pr Re pc u 
P 


St = 


(16) 
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where the total temperatures 


1 2c 


and 


(22) 


2 

2 2c 


The final six differential equations describing the flow characteristics 
derived for the regions of interest are presented in non-dimensional form 
as follows. 

Impingement region, 


1-M du IT d r T p V 

;ir ^ <57 - 

1 


1-^; “1 




”t ymJ/2 s 

'*V"i ' 


r '“1 ^ k 


d^T. 


Y-1..2 dr D.u.c T. ,2 
I i 1 1 p t^ dr 


(23) 


dT. 


T dr 


ird.r T 

(f- - DSt + ( 


Pl"l 


) [7 


dA. 
A^ dr 


(y-i)m: 


1 




. Y-1.. 2 P.u.c .7 

1+2^1 1 P ^ 


(24) 
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M_ 


«2 2 
n r 

yRt. 


^2 2 
n r 

yRT^ 


(1 


Y-1 -2, 

^ V“ 2 > 


(29d) 


where M = Mach number. 

The measured total temperatures and the static pressures at the passage 
hub of both the impingement and supply channels are prescribed as initial 
conditions for Integration of these equations. The velocities at the 
passage hub of the regions are prescribed as u^ = 0 and U 2 = &/p 2^2 where 
m is the total mass flow rate entering the supply region. 
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IV. DATA ANALYSIS 

The solutions of the differential equations derived In III enable one 
to obtain the local flow distribution, which further gives the circumferen- 
tially averaged flux to the coolant along the radius. In this section, a 
method of determining the local flow distribution and heat transfer 
distribution is outlined. 

From the measured heat fluxes Kreatsoulas and Preiser calculated, the 
local heat transfer coefficient h of the leading edge surface using the 
definition 

h = q/(T - T^) (30) 

c c 

where T^ is the measured inlet coolant temperature, and q and T are the 
measured local heat flux and skin temperature, respectively. The value of 
h^, calculated using Eq. (30), is based on coolant inlet temperature to 
the blade, so it includes a number of complex effects, such as heating of 
the fluid in the supply passage and variations in jet Reynolds number due 
to differential pressure drops in the supply and impingement passages. A 
more useful correlation of the local heat transfer distribution may be 
obtained by basing h on the local fluid temperature Tf in the jets, and 
by correlating in terms of the local jet Reynolds number. Thus, 

h^ = q/(T - T^) (31) 

For the channel flow analysis, the heat transfer coefficient is averaged 
around the circumference of the passage along the radius, and Eq. (30) 


becomes. 
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Chupp, et al. (Eq. 34) which Is expressed as a function of geometric 
parameters and the jet Reynolds number. 


XT « r. 0.7f^d) 0.5 

Nu = 0.63 Re^ 
avg d 


0.6 


exp[-1.27l-J|- 


Sj fdl 0.5M1 1.2 


( 34 ) 


Nu is an arithmetic average of the Nusselt number for strips located in the 
leading edge region. Re^j is the jet Reynolds number which is defined as 
Re = VjP 2 d/ii where d is the hole diameter. The friction factor is also 
recalculated as a function of local Reynolds number. 


f 


0.0791 


Re 


0.25 


(35) 


An approximate flow distribution is determined from the calculated 
average Nusselt number and friction factor based on Eqs. (34) and (35), 
respectively, with the prescribed initial conditions. 

The calculated and measured values of the pressure ratio between 
Impingement and supply regions at the inlet are compared. Then an iteration 
method is applied to satisfy the continuity and zero mass flow at the tip 
of the supply region. 

In the next and final iteration, the Stanton number defined in Eq. (16) 
is calculated from the average heat transfer, calculated by using Eq. (32). 
Subsequently, the new local heat transfer averaged in the theta direction, 
which is based on the temperature difference between the local skin 
temperature and fluid temperature, is calculated as defined in Eq. (33). 

This analysis is applied to all the test data obtained by Kreatsoulas 
and Preiser. The new calculated heat transfer distribution is then compared 
with the correlations obtained by other investigators to ensure that the 
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V. RESULTS AND DISCUSSION 

The Independent parameters of primary interest in the present study 
are Reynolds number, rotational speed, heat input, and stagger angle. 

Tables 2 and 3 summarize the series of parametric studies performed by 
Kreatsoulas and Preiser which are considered in this analysis. Figures 6A 
to 34C are the graphical representation of the calculated local flow 
distribution for all the cases. In Figures 6A to 22C, the velocity, 
pressure, and temperature of the fluid in both the impingement and supply 
channels are plotted for each teat conducted at -30 degree stagger angle. 

The results of the local flow distribution calculated for 0 degree stagger 
angle are similarly Illustrated in Figs. 23A to 34C. On the temperature vs 
radius graphs, the temperatures measured by themocouples at the passage hub 
and tip are also plotted. 

In the impingement passage, the initial condition was applied that 
should be equal to the value given by the thermocouple at the base of the 
impingement passage, and in the region below the first jet, conduction in 
the radial direction was Included in the energy equation. This is why the 
temperature shows an initial drop before rising due to convective heating. 

The outer thermocouple gave fluid temperatures quite different from 
those computed, and quite inconsistent with the measured heat transfer 
values. We infer that such measurements are unreliable in the complex, 
thermally stratified flow found in the impingement passage. 

Figures 6D to 22D and 23D to 34D present the Nusselt numbers averaged 
around the circumference of the passage along the radius for the 30 and 0 
degree stagger angle tests. Rotation Introduces important radial variations 
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is subject to minimal cross flow and buoyancy effects, it is dominated by 
the tangential velocity effect and deflects toward the hub. The third jet, 
under the influence of Increasing buoyancy and cross flow, is deflected 
toward the tip and hence creates a large thermal gradient. 

The effect of Reynolds number on heat transfer can be examined by 
comparing Figs. 13D, 15D, and 17D. In all three cases, the tests were 
conducted at medium rotational speed and high temperature ratio at 30 
degree stagger angle. As the average jet Reynolds number varies from 17000 
to 74000, the heat transfer rate Increases about three times, the average 
values ranging from 100 to 280. Similar trends are shown for the low and 
high speed cases but, at the high rotational speed, the thermal gradients 
between the jet holes along the passage are more pronounced. For the 0 
stagger cases. Figs. 26D and 3 ID show the increase in heat transfer rate by 
a factor of 2 as the Reynolds number increases from low to medium. In all 
the cases, the temperature ratio between the average skin temperature and 
the coolant temperature seems to have almost no effect on Nusselt number as 
shown in Figs. 21 D and 22D. 

The stagger angle effect on impingement-cooled rotating blades is 
examined by comparing the data obtained by Kreatsoulas and Preiser, which 
are conducted at 30 and 0 stagger angles, respectively, holding all other 
nominal conditions the same. Figures 220 and 34D show that the heat 
transfer rate are lower at zero degrees by approximately 30%, and the high 
thermal gradient effect due to the dominant tangential velocity effect 
between the second and third jets is not seen at this angle. 

Also Indicated in Figs. 6D to 34D are the comparable Nusselt number 
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VI. CONCLUSION 

A channel flow model of flow in the Impingement cooling passage of a 
rotating turbine blade leading edge is used to correct measured heat 
transfer coefficents to local fluid conditions. The flows in the supply 
and Impingement regions are both modeled as steady state channel flows, 
using friction and heat transfer coefficents. The local flow properties 
are obtained by solving the momentum, mass and energy equations with the 
prescribed initial conditions, and using the local heat transfer rates 
calculated from the data of Kreatsoulas and Preiser. The local Nusselt 
number, based on calculated local fluid temperature, is re-evaluated by 
circumferentially averaging the heat transfer rate. When this result is 
compared with the Nusselt number, based on the inlet coolant temperature, 
the maximum difference between them is determined to be about 17 percent. 
By comparing the calculated Nusselt numbers with the correlations obtained 
by other investigators, it is concluded that the experimental and modeling 
formulations are satisfactory. The modeling and the analysis Illustrated 
here, therefore, offer a means for more accurate prediction of impingement 
cooling heat transfer rates. 
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TABLE 1 : TEST SECTION GEOMETRY 


Span 

S 

101.6 

mm 

4.000 

in 

Hub radius 

Rh 

406.4 

mm 

16.000 

in 

Tip radius 

Rt 

508.0 

mm 

20.000 

in 

Leading edge diameter 


12.7 

mm 

0.500 

in 

Stagger angle (from axial) 

Y 

0.523 

rad 

30.0 1 

deg 

Impingement hole diameter 

d 

2.1 

mm 

0.081 

in 

Impingement hole pitch 

P 

6.1 

mm 

0.240 

in 

Wall to jet distance 

z 

4.1 

mm 

0.162 

in 

Impingement Insert span 

Si 

76.2 

mm 

3.000 

in 
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TABLE 2: -30 degree stagger angle tests 


TEST If 

ROTATIONAL SPEED 

RE If 

TW/TC 

43 

low 

low 

low 

. 42 

low 

low 

high 

39 

low 

medium 

low 

44 

low 

medium 

high 

40 

low 

high 

low 

41 

low 

high 

high 

57 

medium 

low 

low 

56 

medium 

low 

high 

55 

medium 

medium 

low 

. 54 

medium 

medium 

high 

59 

medium 

high 

low 

58 

medium 

high 

high 

63 

high 

low 

low 

65 

high 

medium 

low 

64 

high 

medium 

high 

66 

high 

high 

low 

67 

high 

high 

high 


TABLE 3: 0 degree stagger angle tests 


TEST If 

ROTATIONAL SPEED 

RE If 

TW/TC 

113 

low 

low 

high 

114 

low 

medium 

high 

115 

low 

medium 

high 

123 

medium 

low 

high 

124 

medium 

low 

hifeh 

122 

medium 

medium 

low 

119 

mediiim 

med ium 

high 

120 

medium 

medium 

high 

121 

medium 

medium 

high 

128 

high 

med ium 

low 

126 

high 

medium 

high 

127 

high 

high 

high 
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?ig. 2: Blade Model Cross Section (Side View) 

1. T * Thermocouple 

2 . P <■ Pressiire Tap 


* 

Thermocouple and pressure tap locations 

1. 0.3905 m 

2. 0.4747 m 

3. 0.4008 m 

4. 0.4921 m 


* measured from the center of the shaft 
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Fig 5: Tes-t Section Geomatrv ( Top View) 

1. Supply Plenum 

2 . Jet Hole 

3. Impingement Space 

4. Impingement Insert 

5 . P.ubber Seal 

6 . Cover 

7. Resistive Wall 
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